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ABSTRACT 

As genetic diversity is a valuable tool for plant breeding and screening promising genotypes. The increase in 
population and the subsequent rise in the demand for agricultural produce are expected to be greater in regions where 
production is already insufficient, in particular in South Africa and South Asia. The necessary increase in agricultural 
production represents a huge challenge to local farming systems and must come mainly from increased yield per unit 
area. To meet this requirement various crop improvement programmers have been initiated around the globe. Under 
any crop improvement programme a sample of promising genotypes are performance tested every year. To obtain good 
samples of genotypes many researchers studying the genotypic variation between/ within genotypes. In this addition to 
this chain an attempt is to make to evaluate the loss of information during the study of these genotypes using the 
conditional reference prior. It is observed that that an integrated likelihood approach for studying genotypic variation 
obtained by using a conditional reference prior is equal to the marginal likelihood approach for parameter obtained by 
the non-centrality parameter. 
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INTRODUCTION 

As increase in population and the subsequent rise in the demand for agricultural produce are expected to 
be greater in regions where production is already insufficient, in particular in South Africa and South Asia. The 
necessary increase in agricultural production represents a huge challenge to local farming systems and must come 
mainly from increased yield per unit area. To meet this requirement various crop improvement programmers have 
been initiated around the globe. Under any crop improvement programme a sample of promising genotypes are 
performance tested every year.As genetic diversity is a valuable tool for plant breeding and screening promising 
genotypes. Many researchers focused their attention to study genotypic variation over the last five decades years 
and explode the various concepts relating with genotypes. But the concept related with Bayesian inference for 
genotypic variation yet not is focused. 

Bernardo (1979) initiated the reference prior approach to development of non-informative priors. Because 
of their importance scores of applications and reference prior methods have been progressively defined and refined 
over the subsequent years. Many research articles recorded the evolutions in the methods including Berger and 
Bernardo (1989, 1992a, 1992b), Berger, Bernardo and Mendoa (1989) and Ye (1990). Due to extensive range of 
applicability, an attempt is made to study the concept of conditional reference prior of non-central parameter in 
genotypic variation. It is observed that, in case of marginal experiments, which provides the quantity F for 
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genotypic variation, discarding all the information about the single mean genotypes. It is then apparent that there is a loss 
of information. Here an attempt is made to clarify such loss of information. It has also been shown that the conditional 
reference prior for the parameter C, obtain from the non-central parameter under the analysis of genotypic variation has been 
same when it is obtain through marginal likelihood approach and integrated likelihood approach. 

Marginal Likelihood Approach 

Let yijbe the observation relative to the j th replication of i th genotype. The terms p , p 2 , p k denote unknown 
parameters and Gij represents the random error component, which we suppose to be independent and identically distributed 
(IID) according to the low N(0, a 2 ) with unknown. Therefore 

Yij = + £ ij (i = 1.2, ... K, j = 1,2, ... r s ) 

and Null hypothesis H 0 : p = p 2 =, ... , = p k = p (Say) 

Against Alternative hypothesis p ; ^ p 

based on the test statistic 



£i r i(yi~y) 2 / u 

Si2j(yij-yi) 2 / Ul 
/ u 2 



= k — 1, u 2 = N — k 



Whereyij = y ij; y = Yij and N = £i L 

It is well-known that under the null hypothesis H 0 , the F ca i follows a F u u ^ distribution. Under the alternative 
hypothesis, the follows non-central distribution, whose non-centrality parameter isli.e. 

*■ = ^ZiWj(Pi -p) 2 

Wherejl = Ei w iFi and w, = ^ 

The above testing procedure is a clear and classical example of how the problem of eliminating nuisance 
parameters is handled when the source of variation are genotypes. Although the model has k + 1 parameters, a scalar test 
statistic is contracted to compare different values of the scalar non-centrality parameter X. In a certain sense, the classical 
test acts as if we would have observed only the marginal experiment, which provides the quantity F for genotypic variation, 
discarding all the information about the single mean genotypes (yi,y 2 , ...,yk). 

The great simplicity of the standard one way classification on transformation of multi-parametric problem into an 
equivalent scalar one, where the only parameter involved is the non-centrality parameter, therefore we can rewrite our 
hypothesis as 

H 0 : 1=0 against H-,^ : X > 0 

Then, frequency-based solutions to the problem consist only in the computation of sample realization F ca i of the F 
statistic. Bertolino et al. (1990) propose to consider as the actual result of a marginal experiment and then to use the 
sampling distribution of the F statistic as a marginal likelihood. They also suggest that it is more convenient to work with 
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x 

the parameter^ = - instead of when X — 0(underH o ) the distribution reduces to the well-known Fisher-Snedecor law. 



Therefore the marginal likelihood function associated toF cal is given by 

L rn oc y “ n- r ( 2 +i ) ( UlFcal V 

fy h]~o Pj v 2 7 rf^+il Vu 2 +uiF ra J 



where Pj(Q = exp ^ z ^ z 



If M(a, b; z) is the Kummer’s confluent Hypergeo metric function, then L(Q is given by 

l (0 « (f) (^), (f )' 1 

oc exp 

Where (a)j = is the Pochhammer’ssymbol. 



r-Nq 


|m| 


fu 1 +u 2 u 2 n£ , 


f u lFcal. 


l 2 J 


[ 2 ’ 2 ’ 2 ’ ' 


^u 2 +u 1( F ca ,yj 



By using the asymptotic approximation for L(Q the non informative prior for the parameter C, is [Solari(2002)] 

71(0 a i 

Integrated Likelihood Approach 

Let l p and 0 p denote the p X lvector of l’s and 0’s resp., then the model can be rewrite as 
y = Ap + e 

Where the A matrix is order of n X k and given by 

In 

A = 

OlM k 0[ 

And e is N(0, <r 2 I n ). Then, the parameter ^can be rewrite as 
S = £ (M- - Plk) 1 ^^ - Plk) = £ p'D'RDp 
whereR = idiag(N 1 ; N 2 N k ) 



Ni 


ON! ’ 


o Nl 


N k 


°N k • 


■ °N k . 



D 



N - Ni -N 2 -N k 

-N x N-N 2 -N k 



-N x N 2 N-N k J 

Also we can express the parameter of interest (Q as 



^ = j^2]i Ni Ri(Ri-R) 



Or 
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? = Hlk) = 

Since the k-dimensional vector (P - plk) lies in a (k-1) dimensional subspace of R k , so we need to transform 
(p — ]31| ( ) into a (k-l)-dimensional vector© = (0 1( 0 2 , ... , 0 k _i)'. As rank (RD) =rank (D) =k — 1, there exists a (k — 1)X 
k matrix Psuch thatPP' = RD. Then the Jacobian matrix J x of the transformation from (p, a 2 ) to (0, a 2 , p) is 

PR 1 0 k _i 

o; 1 

Ik o 

Where p = Rtp'G + pl k 

Now we can reparameterizeO in the new parameter 
R = 9'0 = NjCPi - p) 2 
and 




ip = ( 0,n ) k ~ 3 X (0,27r)[Berger et al. (1998)]. Then 



d(0,a 2 ,H) j 



r f' 2 b' 0 0 
Ok-1 1 0 

O'k-1 0 1 

ji 1/2 c 0 k2 0 fc2 . 



Where the ( k — l)-dimensional vector b and the (k — 2) x (k — 1) matrix xc are given by 



b = ri 1 ' 2 



dd 

dr] 



- 1/2 ^0 
C — Tj A/z — 

1 dxh 

Since C = “ 7 - the Jacobian matrix / 3 of the transformation from (r], a 2 ,Ji, ip) to ((, a 2 , Ji, ip) is given 



_ d{T],a 2 ,n,\li) 1 

j3 ld(<;,<T 2 ,ji,w\ 



a 2 0 



0 0' t 



*2 

c 1 0 o'k-2 

0 0 1 0' fc _2 
l 0(£-2 Ofc - 2 0 fc _2 I k —2 J 



Here we defined a new set of parameters ((, a 2 ,p, ip) in which ( is the parameters of interest and f = (a 2 , /i, ip) is 
the nuisance parameters. For this new parameters ({, f) is given by 

Itf.O = 



Where /(p, a 2 ) — 



\a'a 

G* 4 


o k 


= N 




O' 


N 


o' 1 


u k 


2(- 2 ) 2 J 




U o 

2( ff 2) . 
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/ 


rft ft 


ft'ft 


\ 


/ 


f 


O ' 2 


1 ,\ 


/«", f) = Ndiag 


ft'ft 


<6'6+| 




V 


O ’ 2 


(a 2 ) 2 


J 



= Ndiag 



J_ 1 

4£ 4 <7 2 

1 <+2 

4<r 2 



4(- 2 ) 



-4,Ccc' 



Because b’b = ■ 



= Ndiag 



l ll l 12 

hi h2 



• hi I ^44 j(Say) 



Now follow the reference prior algorithm [Berger and Bernardo 1992], the conditional reference priors for ((, f) 
are given by 



-2 



k I *44 1 1/2 K 

TT 1 

7r r/^' 0 ' 2 ) K I133I 2 K 1 
0(1 I ‘22 1 1/2 

W(0 K \h2\~ 1/2 \hih2 - h2hi\ 1/2 * 

If we constructed the reverse reference prior, where the reverse ordering of the parameter (ip,jl,<j 2 ,Q or (<f*, () 
the prior distribution of the parameter of interest £ conditionally on the nuisance parameter <f * is given by 



/((X) = Ma fl (/ 44) i3 3 |“ l } | 2 ]) 



/ 


<+2 1 


\ 


I , 1 


4( ff 2) 2 40-2 

1 1 


\ 


Ucc,^, 




V 


. 4ct 2 4^ . 


/ 



Therefore conditional reference priors are 

(%0 « iini i/2 « ^ 

( ff 2 /l|/’p) 0C 1 * 22 ^ 0C ~2 

^ (%) ^ I * 33 I V2 K 1 

n*(v|/) « n 44 i 1/2 « nlr 1 3 (sim|j i ) k - i - 2 
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CONCLUSIONS 

This shows that an integrated likelihood approach for studying genotypic variation obtained by using a conditional 
reference prior is equal to the marginal likelihood approach for parameter It is also found that when C, is the real quantity 
of interest the marginal likelihood is a correct report of the information available. 
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